DETERMINATION AND USE OF SPECTRAL PEAK INFORMATION AND 



INCREMENTAL INFORMATION IN PATTERN RECOGNITION 



Cross Reference to Related Applications 

5 This application claims the benefit of United States Provisional 

Application Number 60/219,387, filed July 19, 2000, and United States Provisional 
Application Number 60/198,565, filed April 20, 2000. 

Field of the Invention 

10 The present invention relates to pattern recognition techniques and, more 

particularly, relates to the determination and use of spectral peak information and/or 
similar incremental information in pattern recognition. 

Background of the Invention 

15 In pattern recognition, pattern data is usually quantified by converting the 

data into feature vectors. For instance, in speech recognition, speech data is converted to 
feature vectors that commonly have 39 or 40 elements. The feature vectors are 
subsequently used to analyze the pattern data to determine patterns in the data. Generally, 
during training, a number of classes are developed. During real-time processing, new 

20 feature vectors, created from pattern data, can be assigned to particular classes and 
processed accordingly. 

Researchers are constantly trying to decrease the error rate of pattern 
recognition. One way to do this is to modify the number and type of features in a feature 
vector. Some of these changes have improved pattern recognition, and others have not. 

25 Unfortunately, the only way to currently determine if adding features to or 

changing features in a feature vector affects the error rate of a pattern recognition system 
is to use unmodified feature vectors, determine the error rate, modify the feature vectors, 
and determine the new error rate. If the new error rate is better than the original error rate, 
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the additional or changed features have helped pattern recognition. This can be a time 
consuming and laborious process. 

A problem associated more specifically with cepstral features in speech 
applications is that these features can get corrupted by wide-band noise. Thus, the noise 
5 immunity of cepstral features, in feature vectors, are less than ideal. 

For speech applications, one feature that speech systems can analyze is a 
"formant" feature. Voiced sounds have a particular formant structure when viewed in the 
frequency domain. This formant structure is basically a spectral envelope that overlies an 
underlying speech amplitude curve, and it usually has three "humps" that decrease with 

10 increasing frequency. Conversely, unvoiced sounds have a fairly random structure when 
viewed in the frequency domain. Some speech processing systems try to determine 
representative formant features, which can include determining multiple peaks of the 
formant structure. Multiple peak selection can be fairly complex. Most speech processing 
systems also try to determine formant features even in unvoiced speech, which do not 

15 have formant structures. This can make the formant features very noisy for these 
unvoiced speech regions. 

Consequently, what is needed is a better way of overcoming the problems 
of non-ideal pattern recognition when using feature vectors, lengthy and complex 
determination of whether new or different features improve pattern recognition, noise 

20 resistance of feature vectors, multiple peak selection for formant structures, and noisy 
formant features for unvoiced speech regions. 



Summary of the Invention 

Generally, the present invention determines and uses spectral peak 
information, which preferably augments feature vectors and creates augmented feature 
vectors. The augmented feature vectors decrease errors in pattern recognition, increase 
noise immunity for wide-band noise, and reduce reliance on noisy formant features. 



YOR920000199US3 



-2- 



Illustratively, one way of determining spectral peak information is to split pattern data 
into a number of frequency ranges and determine spectral peak information for each of 
the frequency ranges. This allows single peak selection. All of the spectral peak 
information is then used to augment a feature vector. Another way of determining spectral 
5 peak information is to use an adaptive Infinite Impulse Response filter to provide this 
information. 

Additionally, the present invention can determine and use incremental 
information. The incremental information is relatively easy to calculate and helps to 
determine if additional or changed features are worthwhile. The incremental information 
10 is preferably determined by determining a difference between mutual information 
(between the feature vector and the classes to be disambiguated) for new or changed 
feature vectors and mutual information for old feature vectors. 

A more complete understanding of the present invention, as well as further 
features and advantages of the present invention, will be obtained by reference to the 
15 following detailed description and drawings. 

Brief Description of the Drawings 

FIG. 1 is a block diagram of an exemplary system for processing speech in 
accordance with one embodiment of the present invention; 
20 FIG. 2 is a graph of band pass frequency ranges for three band pass filters 

in accordance with one embodiment of the present invention; 

FIG. 3 is a block diagram of an adaptive Infinite Impulse Response (IIR) 
filter in accordance with one embodiment of the present invention; 

FIG. 4 is a flow chart of a method for determining and using spectral peak 
25 information in pattern recognition, in accordance with one embodiment of the present 
invention; 

FIG. 5 is a flow chart of a method for determining and using incremental 
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information in pattern recognition, in accordance with one embodiment of the present 
invention; and 

FIG. 6 is a block diagram of a system for determining and using spectral 
peak information and incremental information in accordance with one embodiment of the 
5 present invention. 



Detailed Description of Preferred Embodiments 

Referring now to FIG. 1, a speech recognition system 100 is shown that 
determines and uses spectral peak information and incremental information in speech 

10 recognition, in accordance with one embodiment of the present invention. System 100 
comprises a speech utterance preprocessor 110, a number of band pass filters 120, 130 
and 140, each producing a band pass filtered output 125, 135 and 145, respectively, a 
number of adaptive Infinite Impulse Response (OR) filters 150, 160 and 170, a feature 
extractor 180, an augmenting device 185, and additional speech processing elements 190. 

15 The speech utterance preprocessor 110 receives speech and generates 

representative speech waveforms, i.e., speech data. The speech utterance preprocessor 
110 may include, for example, an audio-to-analog transducer (microphone) and an 
analog-to-digital converter that respectively transduce the utterances into an electrical 
signal and then convert the electrical signal into a digital signal representative of the 

20 uttered speech. Further, the preprocessor 1 10 may sample the speech signal and partition 
the signal into overlapping frames so that each frame is discretely processed by the 
remainder of system 100. The output signal of the preprocessor 110 is the sampled 
speech waveforms or speech data that is recorded and provided to band pass filters 120, 
130 and 140 and feature extractor 180. 

25 The band pass filters 120, 130 and 140 are preferably selected to have 

non-overlapping band pass frequency ranges. Frequency ranges for three exemplary band 
pass filters are discussed below in relation to FIG. 2. Any number of band pass filters may 



YOR920000199US3 



-4- 



be used, but normally one to three band pass filters will be used, with each band pass 
filter selected to have a band pass frequency range corresponding to a formant frequency 
range. Common formant frequency ranges are 280-710 Hertz (Hz), 870-2250 Hz and 
2250-2890 Hz. By using these ranges, each adaptive IIR filter (150, 160 and 170) should 
5 track one specific peak. Each band pass filter 120, 130 and 140 produces a band pass 
filtered output 125, 135 and 145, respectively. 

These band pass filtered outputs 125, 135 and 145 each feeds into a 
respective adaptive IIR filter of the set of adaptive IIR filters 150, 160 and 170. Each 
adaptive IIR filter tracks a spectral peak and produces information about the peak. A 

10 preferred adaptive IIR filter is discussed below in reference to FIG. 3. Spectral peak 
information produced by each adaptive IIR filter is preferably spectral peak locations si 
through s z and/or spectral peak energy ei through e z . These are discussed in more detail 
below in reference to FIGS. 3 and 4. 

The feature extractor 180 receives the speech data and, as is known in the 

15 art, extracts spectral features from the signal at regular intervals, for example, about every 
10 milliseconds. The spectral features are in the form of feature vectors that are then 
passed on to augmenting device 185. As is also known in the art, the feature vectors 
representing training data may be used to train acoustic models such as, for example, 
band-quantized Gaussian mixture models, which are then used by the system to decode 

20 speech data received during the course of a real-time application. 

Augmenting device 185 concatenates the spectral peak information, {(si, 
ei) (s 2 , e 2 ) . . . (s z , e z )}, onto a feature vector, {(ft, . . . , f n )}, to create an augmented feature 
vector, {(ft, . . . , f n ) (si, ei) (s 2 , e 2 ) (s z? e z )}. The augmented feature vector is then passed to 
the additional speech processing elements 190. Such speech processing elements (not 

25 shown), as is known in the art, could include a linear discriminant analysis routine, 
hypothesis search routines, an acoustic vocabulary, a language model vocabulary, 
acoustic models, and a language model. 



YOR920000199US3 



-5- 



As discussed in more detail below, spectral peak information helps to 
increase the noise immunity of feature vectors and this particular method of estimating 
them also reduces the effect of noisy formant features. 

While the present invention is primarily illustrated herein in the context of 
5 speech recognition, aspects of present invention are applicable to any pattern recognition 
system that uses feature vectors. For instance, incremental information can be determined 
for any system that uses feature vectors. As is known in the art, a feature vector is a 
vector that contains information about a class. A feature vector can be discrete, 
continuous, single dimensional or multidimensional. Other aspects of the present 

10 invention are applicable to any pattern data that are expected to have a spectral peak or 
have single spectral peaks in certain frequency ranges. For instance, an adaptive IIR filter 
can be used to determine spectral peak information for a single spectral peak or a number 
of adaptive IIR filters can be used to determine spectral peak information for single 
spectral peaks in each of a number of frequency bands. Still other aspects of the present 

15 invention are applicable to any pattern data that are expected to have a spectral peak or 
have single spectral peaks in certain frequency ranges, and where the pattern recognition 
system uses feature vectors. For instance, an adaptive IIR filter can be used to determine 
spectral peak information in a certain frequency band and this spectral peak information 
can augment conventional feature vectors. 

20 Referring now to FIG. 2, a graph of the frequency responses for three 

different band pass filters are shown. These band pass filters approximate the formant 
frequency ranges of 280-710 Hz, 870-2250 Hz and 2250-2890 Hz. For instance, band 
pass response 220 approximates the 280-710 Hz formant frequency range, band pass 
response 230 approximates the 870-2250 Hz formant frequency range, and band pass 

25 response 240 approximates the 2250-2890 Hz formant frequency range. By using these 
ranges, there should only be one spectral peak (corresponding to an associated formant 
peak) in each of the frequency ranges. This allows an adaptive IIR filter to track, over 
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time, the one major spectral peak in its associated band pass frequency range. 

Turning now to FIG. 3, an exemplary adaptive IIR filter 300 is shown. The 
adaptive IIR filter 300 comprises a number of adders 310, 320, 330, 340, and 350, two 
delays 360 and 370, and a number of coefficients k i5 2, and G. Assuming that there are 

5 three band pass filters and denoting the outputs of these three band pass filters 120, 130 
and 140 (see FIG. 1) as y } (t), y 2 (t) and y 3 (t), the adaptive IIR filter 300 identifies the 
spectral peak in the band-limited spectra of yt(t). 

This adaptive filter structure has previously been described in 
Padmanabhan et al., "Resonator-Based Filter-Banks for Frequency Domain 

10 Applications," Institute for Electronic and Electrical Engineers Trans. Circuits and 
Systems, Oct. 1991, the disclosure of which is incorporated by reference herein. 
Generally, the adaptive filter is a multiple notch IIR filter with the notch frequencies 
being directly related to the filter coefficients (each notch frequency depends on exactly 
one filter coefficient). This multiple notch transfer function is obtained by embedding 

15 several digital resonators in a feedback loop, with the notch frequencies corresponding to 
the resonator frequencies. The goal of the adaptive algorithm is to minimize the power of 
the output of the notch filter. For the case where there are N sinusoids in the input, and 
there are N notches in the filter transfer function, the optimal solution is when the notch 
frequencies are equal to the input sinusoidal frequencies. In "Resonator-Based 

20 Filter-Banks for Frequency Domain Applications," an adaptive algorithm was described 
that guarantees convergence under certain conditions, and that has a complexity that is 
linear in N. 

For use with the present invention, it can be assumed that each band pass 
output, y t (t) 9 has a single spectral peak. The adaptive IIR filter 300 of FIG. 3 has a single 
25 notch and can be used to track each of the y l (t). The transfer functions from the input to 
various nodes in the filter structure 300 are given below: 

TJ ( \ _ _*£_ „ l-(2-^)2- 1 +2- 2 
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ttfb,i(Z)- *m - \-(2-fy(\-G)z-*+([-2G)z-2> K 4 " ' 

#W( Z ) ~ ~ \-{2-k}){\-G)z'^{\-2G)z-^ ^ Eq * 3 ^ 

The transfer function H e (z) represents a notch filter, with the notch 
frequency, m h being related to the filter coefficient, k h through the following equation: 
5 ^2sin(f) . (Eq.4) 

The algorithm for adapting the filter coefficient is given by: 

i / ix 7 / \ *e{n)x pSyl (n) . 

+ l) = ki(n) -fl ^(n^n^e ^ 

The term x e (ri)x ps ,i(ri) denotes the pseudo-gradient of the objective 
function (i.e., the coefficient is adapted in a direction opposite the pseudo-gradient), and 

1 0 the term < x pSil (n)x ps j(n) > +e denotes a power normalizing term that modifies the gradient 
direction to point in the Newton direction (i.e., the normalizing term approximates the 
inverse of the Hessian of the objective function). The e term prevents a divide-by-zero 
operation, should the average power term get very small The spectral peak location can 
be inferred from the value of the filter coefficient k h after the coefficient has converged, 

15 using Equation 4. The ju is a rate-determining parameter. As is known in the art, one can 
control the rate at which the filter coefficient changes. The rate-determining parameter 
allows one to control the rate of change of the filter coefficient, and the parameter is 
usually a heuristic parameter. It is determined with sample data, and a variety of 
parameters are tried, with the best parameter selected. The selected rate-determining 

20 parameter is then used during real-time analysis. 

An additional feature of the filter structure 300 is that the transfer function 
Hjbj is the complement of the notch transfer function (H e ) 9 i.e., it represents a band pass 
transfer function with the center frequency corresponding to the resonator frequency, mi. 
Consequently, the power of the signal at x^j represents the power of the input speech 

25 signal at this frequency. 

Thus, the spectral peak information produced by adaptive IDR. filter 300 
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comprises the spectral peak location, Wi (denoted as Si, s 2 and s z in FIG. 1), and the 
energy of the signal at x pj (denoted e u e 2 and e z in FIG. 1). A benefit of this filter 
structure is that this spectral peak information is less susceptible to wide band noise. 
Wide band noise should not change the location of a spectral peak in a formant frequency 
5 range, but could change the level of the peak. Thus, the spectral peak location should be 
relatively immune from wide band noise. The energy at the spectral peak location should 
also be relatively immune from wide band noise because the energy is determined by a 
function that is centered at the spectral peak location and that has a very small amount of 
bandwidth. Therefore, the wide band noise would have to be fairly concentrated at the 

10 spectral peak for the energy to be erroneously high. 

Another benefit of the adaptive HR filter 300 is that the spectral peak 
information should stay on the same trajectory. For instance, if the filter coefficient 
converges to a particular number (indicative of a spectral peak at a particular frequency, 
wi) during a time period with voiced speech segments, and then unvoiced speech 

15 segments appear at the input x h it will take a while for the filter coefficient to change to a 
new spectral peak frequency of the unvoiced speech segments. In fact, it may not change 
at all, because unvoiced speech generally contains few large peaks to which the adaptive 
HR filter 300 would converge. This delayed action tends to smooth the response for both 
the energy and spectral peak locations. In essence, the spectral peak information is less 

20 arbitrary. This is in contrast to normal formant estimates, which will contain arbitrary 
formant features for unvoiced speech segments, as unvoiced speech segments do not have 
a formant structure. 

It should be noted that any spectral peak information could be used with 
the present invention. The examples of the spectral peak information given above are 

25 preferred, but other spectral peak information could be used. For example, formant peak 
estimates could be used instead of a spectral peak information determined by adaptive HR 
filter 300. If formant peak estimates are used, band pass filters and adaptive HR filters 
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would not be needed. The energy of a spectral peak and the spectral peak location are 
preferred over formant peak estimates because the latter can be noisy if the system is 
affected by wide band noise. 

Turning now to FIG. 4, a method 400 is shown for determining and using 

5 spectral peak information in pattern recognition, in accordance with one embodiment of 
the present invention. Method 400 is used during real-time processing of patterns, but 
could be used after pattern data has been accumulated. Method 400 begins in step 410 
when pattern data are retrieved. This pattern data will generally be digital samples of 
analog data. For example, speech data is generally digital samples of an analog waveform 

10 created by a microphone, while handwriting data might be the digital samples of an 
electrical transduction of a pressure plate. 

After the speech data are retrieved, the pattern data are routed to both steps 
420 and step 440. In step 420, the pattern data is split into a number of frequency bands. 
One option at this point is to split the pattern data into three frequency bands, each 

15 frequency band corresponding to a frequency range that should contain a formant peak. 
This is an excellent option for speech data. Another option is to split the pattern data into 
two or less frequency bands, which could also be used on a speech processing system if 
the frequency bands are chosen to coincide with two of the three formant frequency 
ranges. Yet another option could be to split the pattern data into more than three 

20 frequency ranges, which could be useful for visual pattern recognition. To split the 
pattern data, preferably band pass filters are used and designed with the appropriate band 
pass frequency ranges. 

In step 430, the spectral peak information is determined for each frequency 
band. Generally, the spectral peak location and/or the energy of the spectral peak location 

25 are determined in this step for each of the frequency bands that have been band passed. In 
step 435, the filter coefficients for each of the adaptive HR filters are adjusted. If adaptive 
IIR filters are not being used, then, if needed, other coefficients may be adjusted. For the 
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adaptive IIR filters, an algorithm for adjusting the filter coefficients has been discussed 
above in reference to Equation 5. 

Step 440 is performed in parallel with steps 420, 430, and 435. In step 
440, the feature vectors are determined. This generally occurs through the use of a 
5 "feature vector extractor," which is designed to create feature vectors. As an illustration, a 
speech feature vector extractor could determine feature vectors that have 40 features. In 
step 450, the feature vectors and the spectral peak information are combined into 
augmented feature vectors. Preferably, an augmented feature vector would be a 
concatenation of the spectral peak information and a feature vector. This has been 
10 discussed (and shown) in reference to FIG. 1 . 

In step 460, other pattern recognition steps are performed. For speech, 
these steps could include a linear discriminant analysis routine and a hypothesis search 
routines. The end result of method 400 should be a feature vector that has been associated 
with a class and that is used to determine a pattern. The pattern could be a phone of 
15 speech that is then converted to text or a pattern equivalent of handwritten symbol that is 
then converted to a textual symbol. 

It should be noted that these steps will usually repeat indefinitely. For 
instance, in a speech processing system, new speech data will constantly be input into the 
system. This speech data will be converted to text (or compressed) and steps 410 through 
20 460 will be continuously repeated until the speech processing system is turned off. 

By using the present invention, the feature vectors should be more immune 
to wide band noise and should have elements that better track spectral peaks for data with 
spectral peaks. 

Turning now to FIG. 5, a method 500 is shown for determining and using 
25 incremental information. Method 500 is performed whenever it is desired that a new 
feature be added to a feature vector or a changed feature replace an old feature in a feature 
vector. Method 500 begins when classes are determined (step 510). In many situations, 
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such as for speech processing, class determination can be quite complex and involved. 
For the purposes of method 500, it can be assumed that classes can be properly 
determined. For speech processing systems, the classes that are determined are phonetic 
classes. 

5 In step 520, the current features are determined. In step 530 the "old" 

mutual information is determined. The usefulness of a feature may be measured by the 
amount of information it provides in discriminating between classes. This can be 
quantified by the mutual information between the feature vector and the class. Let c 
denote the class and Z denote the feature vector. The mutual information between Z and c 

10 is defined by: 

/(Z; c) =Zp(c) \ z p(Z/c) log [*jj$]dZ. (Eq. 6) 

Though Equation 6 cannot be expressed in closed form, by vector 

quantizing Z into an appropriate number of codewords and approximating the integral 

with a summation, it may be written as: 
15 /(Z;c) =Zp(c) T,p(Zj/c) log[^f ] (Eq. 7) 

Equation 7 is an explicit relationship between a feature vector and classes. 
It would be helpful to measure the amount of information available by augmenting the 
current feature vector with new features. To do this, the new features are determined (step 
540). If the new feature is written as z, a new feature vector, Z, would be created by 
20 concatenating Z with z (step 550). The new mutual information between the augmented 
vector Z and c can be computed from the following: 

/(Z;c) ^p(c) Xp(Zj/c) logf^]. (Eq. 8) 

Equation 8 is performed in step 560 of method 500. The amount of 
incremental information added by z can be computed through the following: 
25 Slz = I(Z;c)-I(Z;c). (Eq. 9) 

Equation 9 is performed in step 570 of method 500. Determining the 
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mutual information in the prior equations in method 500 is relatively simple and not very 
time consuming. As such, the incremental information determined in Equation 9 provides 
an indication of the benefit of a feature, yet costs little computationally to determine. 
Equation 9 can be used to determine if the mutual information of a new feature vector 
5 meets a predetermined increase in mutual information over the "old" mutual information. 
As an example, if a new feature adds 0.1 bits of incremental information, this could be 
significant if the old mutual information was 1.62 bits, as this would be a 6 percent 
increase in mutual information. 

To use the incremental information from method 500, once one feature is 
10 analyzed using method 500, another new or changed feature can be analyzed by 
re-performing steps 540 through 570 with the new/changed feature. By following this 
process, the features that add the most amount of incremental information can be used to 
determine new baseline feature vectors that can be used in step 520 and compared to 
subsequent, even newer feature vectors in steps 530 through 570. 
15 Although the previous discussion centered on one mutual information 

equation, any technique for determining the mutual information between a feature and a 
class may be used. For instance, in Padmanabhan, "Use of Spectral Peak Information in 
Speech Recognition," National Institute of Standards and Technology (NIST), Speech 
Transcription Workshop, May 16-19, 2000, the disclosure of which is hereby 
20 incorporated by reference, the following equation is developed for mutual information: 
J(Z;c) = log US^ r |-SPdog \A2, e A T \, (Eq. 10) 

which turns out to be almost exactly the object function for the heteroschedastic 
discriminant analysis as described in Saon et al., "Maximum Likelihood Discriminant 
Feature Spaces," International Conference on Acoustics, Speech and Signal Processing 
25 (ICASSP), 2000, the disclosure of which is also incorporated herein by reference. 

Turning now to FIG. 6, an exemplary system 600 is shown that could be 
used to perform the methods and apparatus of embodiments of the present invention. 
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System 600 comprises a compact disk 605 and a computer system 610 that itself 
comprises processor 620 and memory 625. Memory 625 comprises some or all of system 
100, method 400 and/or method 500. As such, memory 625 will configure the processor 
620 to implement the methods, steps, and functions disclosed herein. The memory 625 

5 could be distributed or local and the processor 620 could be distributed or singular. The 
memory 625 could be implemented as an electrical, magnetic or optical memory, or any 
combination of these or other types of storage devices. Moreover, the term "memory" 
should be construed broadly enough to encompass any information able to be read from 
or written to an address in the addressable space accessed by processor 620. With this 

10 definition, information on a network is still within memory 625 of system 600 because 
the processor 620 can retrieve the information from the network. It should be noted that 
each distributed processor that makes up processor 620 will generally contain its own 
addressable memory space. 

It should also be noted that computer system 610 could be an 

15 application-specific integrated circuit that performs some or all of the steps and functions 
of system 100, method 400 and/or method 500. 

As is known in the art, the methods and apparatus discussed herein may be 
distributed as an article of manufacture (such as compact disk 605) that itself comprises a 
computer readable medium having computer readable code means embodied thereon. The 

20 computer readable program code means is operable, in conjunction with a computer 
system, to carry out all or some of the steps to perform the methods or create the 
apparatuses discussed herein. The computer readable medium may be a recordable 
medium (e.g., floppy disks, hard drives, compact disks, or memory cards) or may be a 
transmission medium (e.g., a network comprising fiber-optics, the world-wide web, 

25 cables, or a wireless channel using time-division multiple access, code-division multiple 
access, or other radio-frequency channel). Any medium known or developed that can 
store information suitable for use with a computer system may be used. The 
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computer-readable code means is any mechanism for allowing a computer to read 
instructions and data, such as magnetic variations on a magnetic medium or height 
variations on the surface of compact disk 605. 

In Padmanabhan, "Use of Spectral Peak Information in Speech 
5 Recognition," it is shown spectral peak information increases the mutual information (and 
consequently the incremental information) and can decrease word error rate in speech 
recognition. 

Thus, what has been shown are systems and methods for determining 
spectral peak information of pattern data through adaptive IIR filters, for augmenting 

10 feature vectors with features that are more robust, and for determining if a feature is a 
worthwhile addition to already existing feature vectors. 

It is to be understood that the embodiments and variations shown and 
described herein are merely illustrative of the principles of this invention and that various 
modifications may be implemented by those skilled in the art without departing from the 

15 scope and spirit of the invention. For instance, other spectral peak information, such as 
formant estimates, may be used and a different spectral peak tracking device could be 
used in place of the adaptive IIR filters of embodiments of the present invention. 
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